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account for the well known instability of ordinary pipe flow. 

This report describes a more general theory which shows promise of 
overcoming the above limitations. The new theory involves a number of 
innovations in the formulation of some of the basic equations, in the fori 
lation of certain boundary conditions, and in the formulation of the 
criterion of stability itself. The analytical work is essentially comple: 
and preliminary calculations, while still of limited scope, appear to 
support the theory and are therefore quite encouraging. 

Nevertheless, before the new theory can be considered as definitivel 
confirmed, a considerable amount of systematic computer calculation will 
still be necessary. Research support is now being sought for this purpos 
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SUMMARY 



This report provides a detailed technical outline and evaluation of 
recent research performed at the Naval Postgraduate School on the hydro- 
dynamic stability of plane Poiseuille flow and of pipe Poiseuille flow. 

These two cases are of fundamental importance in connection with the 
theory of hydrodynamic stability. Each case involves significant discrep- 
ancies between the predictions of conventional theory and the results of 
actual experimental observations. In particular, the conventional theory 
fails completely to account for the well known instability of ordinary 
pipe flow. 

This report describes a more general theory which shows promise of 
overcoming the above limitations. The new theory involves a number of 
innovations in the formulation of some of the basic equations, in the formu- 
lation of certain boundary conditions, and in the formulation of the 
criterion of stability itself. The analytical work is essentially complete 
and preliminary calculations, while still of limited scope, appear to 
support the theory and are therefore quite encouraging. 

Nevertheless, before the new theory can be considered as definitively 
confirmed, a considerable amount of systematic computer calculation will 
still be necessary. Research support is now being sought for this purpose. 
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1. Introduction 



The basic research discussed in this report deals with the classical 
problem of the hydrodynamic stability of parallel flows. Two fundamental 
cases, both of considerable historical and practical importance, are of 
particular interest, namely, the plane Poiseuille flow between infinite 
parallel plates and the flow through a uniform pipe. Despite the consid- 
erable amount of research that has been done over the years on these two 
fundamental cases, the theory available up to now does not account ade- 
quately for the observed experimental facts. Indeed, for the case of pipe 
flow, the classical theory fails to predict the onset of hydrodynamic 
instability at any Reynolds number whatever! In this respect the existing 
theory is clearly inadequate. 

A critical review of the situation by the present investigator has dis 
closed that the above limitations and defects of the classical theory can 
largely be overcome by generalizing the analysis to include certain three 
dimensional effects that are neglected in the classical approach. The re- 
sulting calculations have succeeded for the first time in revealing instab i 
lities in pipe flow for certain types of three dimensional perturbations; 
these preliminary results appear to be consistent, at least qualitatively, 
with the known experimental facts about pipe flow. The corresponding 
calculations for plane Poiseuille flow reveal that under certain conditions 
instability can occur at a Reynolds number significantly lower than that 
predicted by the classical theory. The new analysis therefore promises to 
advance our basic understanding of the important phenomena of hydrodynamic 
instability in a significant way. 
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The exploration of the quantitative implications of this generalized 
three dimensional theory entails a heavy burden of computer calculations. 
Research support is needed for the purpose of carrying out systematic 
calculations of this kind. Support is also needed for the purpose of 
documenting and publishing in the technical literature the analytical 
basis of the theory and such preliminary numerical results as have already 
been obtained. 

More specifically, it is proposed to carry out the above documentation 
and publication, as well as a significant amount of the needed additional 
calculations during the first year of the proposed research. If warranted 
by the results accomplished during the first year, additional support for 
follow on research will probably be sought also for another year or two 
thereafter, to consolidate the progress made. 
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2. Analytical Approach 

The essentials of the analytical approach involved can best be ex- 
plained by first considering the specific case of plane Poiseuille flow. 
All equations are nondimensionalized by utilizing as dimensional 
reference parameters the density p of the fluid, the semi-height a 
between the walls and the mean volumetric velocity U of the fluid, 

Cartesian coordinates x, y, z are used with x in the direction of 
the mean flow and y normal to the walls. Dimensionless coordinate y 
varies from -1 at the lower wall to +1 at the upper wall. Symbols 



i, j, k denote unit vectors in the x, y, z directions, respectively. 

It is well known that the dimensionless velocity U of the mean flow 
along the x axis is given by the simple expression 



The hydrodynamic stability of this flow field can be analyzed in the 
usual way by superimposing upon the above mean flow a suitable perturbation 
velocity v of small amplitude. It is of course necessary that v satisfy 
the continuity equation for incompressible flow, namely. 



A convenient way to ensure that Eq. (2.2) is satisfied is to express 

-> 

v in terms of a corresponding vector potential function, let us call it 
-> 

W , by a relation of the form 



It is customary in linear stability analysis to represent quantities 

-T -> 

like v and W by complex variables with the understanding that only 



U = | (1 - y 2 ) 



( 2 . 1 ) 




( 2 . 2 ) 




(2.3) 
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the real part of the variable shall be identified with the actual physical 
effect in question. This convention simplifies the analytical form of 
the equations involved. 

For a given perturbation vector field v , the corresponding vector 
potential function W is not uniquely determined but may be expressed in 
a number of alternative forms all of which are ultimately equivalent. In 
the present case, however, it proves advantageous to represent the complex 
vector potential function W as follows 

w = [i(o) + j G(y) + k H(y) ] e (aX + 02 + yt) (2.4) 

where the complex constants a, B, y may be further expanded to the form 

a = a + i a_ 

K. I 

6 ■ S R + 1 6j (2.5) 

Y = T R + i 

and where i denotes the imaginary unit, that is, 

i = + /-l (2 . 6) 

The quantities G(y) and H(y) in Eq. (2.4) amount to complex stream 
functions. Although a single stream function suffices to represent any 
incompressible plane flow, two stream functions are required to represent 
an arbitrary incompressible three dimensional flow with complete generality, 

Incidentally, this use of two stream functions to represent the pertur- 
bation flow in this problem is an innovation introduced by the present 
investigator. As far as he knows, this method does not appear in the 
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existing technical literature which deals with the present problem. 

Of course, the actual forms of the complex functions G(y) and H(y) 
are initially unknown but are ultimately calculable in the course of the 
final numerical solution. 

The notation used for the exponential coefficients in Eqs. (2.4) and 
(2.5) was adopted for reasons of overall generality, consistency and 
simplicity in the subsequent analysis. Various alternative but equivalent 
notations are also found in the technical literature, but all such alter- 
natives can be readily translated to the above format whenever necessary. 

Parameters ot^ and 6^. in Eq. (2.5) are associated with the 
imaginary unit i , and hence they represent the oscillatory wave numbers 
of the perturbations with respect to space coordinates x and z, re- 
spectively. On the other hand a and B~ do not involve i and hence 
they denote real exponential growth rates (or, if negative, decay rates) 
with respect to space coordinates x and z . Likewise, y represents 
the oscillatory wave number of the perturbation with respect to time t , 
that is, the angular frequency, while y D represents the exponential 
growth or decay rate with respect to time. Incidentally, it is convenient 
to restrict and B to non-negative values and there is no loss of 

generality involved in doing this so long as a and 6 remain free to 

R R 

take on positive or negative values as required. 

Notice that the above quantities and equations are defined initially 

in relation to axes that remain fixed with respect to the walls. Later we 

shall also consider axes that translate at constant velocity relative to 

the fixed walls. Such motion does not affect parameters a , a , 6 , B T 

R X R I 

but can in general change the apparent values of y and y . We denote 

R I 

these changed values with respect to the moving axes by symbols y^ and yj. 



2-3 



In the method of analysis considered here, the numerical values of 



a , a , 0 , 6 are assigned arbitrarily, but the corresponding values of 
K 1 K L 

y and y must be found by calculation. This point will be clarified 

further in the subsequent discussion. 

Of course it is essential that the resultant perturbed flow field 
defined by the foregoing relations satisfy Newton f s second law of motion. 
This in turn requires that the velocity and vorticity components satisfy 
the vorticity transport equation, which is the relation that is found by 
applying the curl operator to the vector equation of motion. This relation 
involves terms which are linear in the perturbation quantities and terms 
which are quadratic in the perturbation quantities. For small perturba- 
tions, however, the quadratic terms may be neglected. In this way the 
customary linearized equation is obtained. A full explanation of these 
details lies outside the scope of the present discussion, but is given by 
Harrison, Ref. (6). Fortunately, the result is well known. One convenient 
version of this result may be summarized concisely as follows 



e 




( 2 . 7 ) 



where the symbols are defined as indicated below. 



U = i | (1-y 2 ) 



Velocity of mean flow field 



ft = VxU 



Vorticity of mean flow field 



v = VxW 



Perturbation velocity 



( 2 . 8 ) 



7 xv 



Perturbation vorticity 



R 



Reynolds number 



e 
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Notice that Eq, (2,7), the vorticity transport relation, introduces 
yet another important parameter into the analysis, namely, the Reynolds 
numb er R 

e 

Of course the solution of Eq t (2,7) must be such as to satisfy the 
boundary conditions of the problem which are simply that 

v = 0 at y = + 1 (2,9) 



By a lengthy but systematic procedure which is fully documented in 
Ref (6), Eq. (2.7) can be expanded in terms of W and hence in terms of the 
functions G(y) and H(y) . Eq. (2,7) is a vector relation and therefore 
implies three corresponding scalar equations. It can be shown, however, 
that only two of these are independent. Consequently, the final problem 
reduces to two differential equations in the two unknown functions G(y) 
and H(y) . 



In this connection, the following three auxiliary functions occur. 
These functions are immediately known as soon as definite values are 
assigned to parameter R^ and to the complex constanst a and 8 . It 
is helpful to calculate these functions first, before proceeding with the 
rest of the solution. Thus 



T(y) - - . f U-y 2 ) 

e 



I'y) = (C * 2 * ^ + T(y) 



( 2 . 10 ) 



J(y) = (a 2 + 8 2 ), T(y) - 3a 



Using symbol D to denote the differential operator — , we find that 
the resulting vorticity transport equation and its associated boundary 
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condition can now be summarized concisely in the form 




+ I(y) 



d 2 h 



+ J (y) H > - y 



D 2 H + 



(a 2 + 



B 2 ) 



H> = 0 



where H = 0 
and DH = 0 



at y = + 1 



( 2 . 11 ) 



The second vorticity transport equation and its associated boundary 
condition can likewise be written 

j|- D 2 G + [T(y) - YJ G 
I e 

= — 2 — ^-~ 2 ~ )|- D 3 H + [T(y) - y] DH - 3ayH[ (2.12) 

(a + £ ) ( e ) 

where G = 0 at y = + 1 

Now consider the situation in which the values of parameters 

a , a , g , 3 t and R are arbitrarily assigned. Then the three auxilary 
K X R X e 

functions defined by Eqs . (2.10) become known functions. Hence the basic 
relations given in Eqs. (2.11) can now be solved numerically, subject to 
the boundary conditions shown. 

Eqs. (2.11) define an eigenvalue problem. Theoretically there exist 
an infinite number of distinct solutions to this equation, each characterize 
by a corresponding value of the complex constant y . We term these roots 

eigenvalues and denote them by symbols 'V'l’ ^2* ^3 Associated with 

each root is a corresponding complex eigenfunction H^(y) , ^(y), H^(y).., . 

Fortunately, only the solution which corresponds to the least stable root 
is required for the stability analysis. The least stable root is defined 
in more detail later. 
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are 



Once the kth root y^ and the kth eigenfunction H^(y) 
known, Eqs. (2.12) may be solved for the corresponding complementary 
eigenfunction G^(y) > subject to the boundary condition shown. 

Notice that Eqs. (2.11) and (2.12) for H and G , respectively, 
are uncoupled and may therefore be solved sequentially. This fact simpli- 
fies the solution. Usually we are interested only in finding the least 
stable root. In that case it is unnecessary to solve Eqs. (2.12). 

By utilizing the method of finite differences, it is possible to 
reduce the foregoing equations to matrix format. The resulting matrix 
represents a standard eigenvalue problem which can be solved numerically 
by familiar and well established techniques. 

The details of this reduction to matrix format are presented in 
Ref. (6) and are not repeated here. 

The real part of each eigenvalue is denoted by y and the imaginary 

R 

part by y^ . For reasons that are more fully explained later, the eigen- 
value that has the algebraically highest value of y is termed the least 

R 

stable root. We denote this highest real value by symbol y . Of course 

RL 

y depends on the particular numerical values assigned to parameters 
RL 

a n , a T , 8 t and R . We may symbolize this relation in the form 

R I R I e 



Y RL Y RL ^R’ a I’ 6 R» 6 I* R e ^ 



(2.13) 



The stability of the flow field is determined by the numerical value 



of y according to the criteria considered in the next two sections. 
RL 
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3. Conventional Criterion of Stability 

Historically, the analysis of stability has concentrated attention 
chiefly on the important but nevertheless special case in which the spatial 
form of the pertubations is taken to be purely sinusoidal. This amounts 
to restricting the analysis to the particular case for which, in the present 
notation, 

°R 2 + ^R 2 = ° (3,1) 

Moreover, any time that Eq. (3.1) applies, the well known theorem of 
Squire, Ref. (1), also applies. When expressed in terms of the present 
nomenclature, Squire f s theorem simply asserts that for specified values of 
parameters and a^. , the most unstable condition will always be the 

one for which 

Sj = 0 (3.2) 

Thus of the five parameters R a , a , B , $ which occur in the 

e K 1 R I 

general case, the constraints expressed by Eqs. (3.1) and (3.2) eliminate 
three. Therefore the conventional analysis finally retains only two para- 
meters, namely, R and a. . 

e I 

In the conventional nomenclature the quantity we have called a^. is 
usually represented by the unsubscripted symbol a ; to avoid ambiguity in 
the present more general discussion, however, it is important for clarity to 
adhere consistently to the subscripted notation as introduced originally in 
connection with Eqs. (2.5). 

Under conditions where Eq. (3.1) applies, that is, for the conventional 
analysis, the stability of the flow field is determined by a very simple 
rule. The flow is stable, or at least neutrally stable, provided that 
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(3.3) 



If the inequality applies, the perturbation amplitude decays expo- 
nentially over time and the flow is stable. If the equality applies, the 
perturbation amplitude neither grows nor decays overtime and the flow is 
neutrally stable. If neither of the above conditions is satisfied, the 
perturbation amplitude grows with time and the flow is unstable. 

Numerical calculations based on the conventional restriction of Eq, 
(3.1), on the application of Squire T s theorem as expressed by Eq. (3.2), 
and on the simple stability rule summarized in Eq. (3.3), permit the flow 
stability to be determined for arbitrarily assigned values of and 

a_£ . The region of instability in the R^ vs a^. plane has the general 
form sketched qualitatively in Fig. 3.1. Notice that there exists a 
critical Reynolds number below which sinusoidal streamwise perturbations 
of all wave numbers a^. are stable. The conventional analysis gives the 
value of the critical Reynolds number, based on channel semi-height and 
volumetric mean velocity as about 3850, Thomas, Ref. (2). Of course, all 
of these conventional results are familiar; they have been known for some 
time . 
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STABILITY OF PLANE POISEUILLE FLOW 



i 









i 



4. Generalized Criterion of Stability 

The simple stability rule summarized in Eq. (3.3) of the previous 
section, although correct for the special case of purely sinusoidal 
oscillations, is not adequate for dealing with the more general situation 
tfhere 



Moreover, under these conditions, Squire’s theorem no longer holds 
either, so that 



Thus of the five parameters R , a , a , B , S T which occur in the 

e K 1 K 1 

general case, none can now be dropped. Compared with the conventional 



i considerable increase in complexity. The additional parameters which 
)ccur in the generalized problem are associated with certain instabilities 
rfiich are not revealed in the conventional analysis. This is the basic 
reason why the generalized analysis gives results in better agreement 
?ith experiment than does the conventional method. 

While the simple criterion of stability expressed by Eq. (3.3) is 
miversally accepted as applying to the conventional case, there is no 
such unanimity among different investigators concerning the appropriate 
stability criterion to apply to the generalized problem. 

The analysis of the previous section has shown that if we arbitrarily 
Impose the restrictions that 



( “r 2 + 6 r 2 > * 0 



(4.1) 



e T t o 



(4.2) 



:ase which involves only the two parameters R^ and , this represents 



a. 



■r = 0 



8 . 



R = 0 



( 4 . 3 ) 
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these conditions are sufficient to ensure neutral stability, Eqs. (4.3) 



are frequently misinterpreted as representing the conditions which are 
necessary to ensure neutral stability. In fact Eqs. (4.3), while un- 
questionably sufficient, are not necessary for this purpose. If they 
were necessary conditions, they would of course rule out of consideration 
altogether the generalized case which is the very subject of the research 
under discussion. Nevertheless, this seems to be a fairly common mis- 
understanding of the situation, a fact which might help explain why the 
fully general three dimensional case has not received the degree of 
research attention that it deserves. 

A useful clue to a more satisfactory formulation of a generalized 
stability criterion is provided by the existance of two special cases in 
which there seems to be general agreement on how the stability should 
properly be defined. Refs. (4) and (5). These are summarized in Eqs. 

(4.4) and (4.5) below. 

The flow field is stable, or at least neutrally stable, 

if « R = o | 

when 1 0 ) 

or if Y r = 0 ) 

when 1. 0 / 

A significant feature of Eqs. (4.4) and (4.5) is that they both refer 
explicitly to parameter a. but not to fi . Thus stability is defined 
in terms of the component of spatial growth rate in the strearawise directio; 
but not in terms of the component in the transverse direction (if any) . 



(4.4) 



(4.5) 
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In analytical work such as considered primarily in this discussion, 
parameter o can be set to zero independently as required by the first 
of Eqs . (4,4) and g and the other parameters can be chosen arbitrarily. 

The value of f or the least stable root can be found by calculation, 

and the resulting stability determined from the second of Eqs, (4.4). 

While this method is favored by some investigators, it fails for the 

important case where a f 0 . 

R 

In a hypothetical experiment involving steady oscillations we would 

have y = 0 in accordance with the first of Eqs. (4.5). The observed 
R 

streamwise growth or decay rate would then fix the value of a for the 

least stable root. The second of Eqs. (4.5) would then fix the stability. 

While this concept can hardly be faulted for the special case to which it 

applies, unfortunately it fails for the important case where y R r 0 . 

One basic reason for the diversity of stability criteria utilized by 

different researchers is that in the general case the apparent exponential 

growth or decay rate of the perturbation with respect to time, as denoted 

by parameter y’ , depends in part on the motion, if any, of the observer’s 
R 

reference frame. To see this consider axes that translate at constant 
velocity parallel to the fixed walls . Let x q denote the translation 
velocity component in the downstream or x direction, and let z Q 
denote the velocity component, if any, in the lateral or z direction. 

Let x T , y, z T denote the coordinates of a point fixed with respect to the 
moving axes, and let x, y, z denote the coordinates of the same point 
with respect to fixed axes. Then 

X = X* + X t 

° (4.6) 

z = z ' + z t 
o 



4-3 



The exponent in Eq. (2.4) can now be written as follows 



X = ax + Sz + yt 



= a(x' + x t) + 3(z' + z t) + yt 
o o 



= ax' + Sz' + (y + ax +3z)t 

o o 



(4.7) 



= ax' + Sz' + y't 



where 



y' = y + ax + $Z' 



(4.8) 



Breaking this into real and imaginary parts gives 



Y R = Y R + “r X o + 8 R 2 o 



(4.9) 



Y i ■ Y I + “i x o + B I Z o 



Thus, with respect to a reference frame which moves with any specified 



values of x and z , the flow field may be said to be stable, or at 
00 



least neutrally stable, provided that for the least stable root 



y 'rl ■ (y el + “A + s rV ^ 



(4.10) 



The least stable root is that root which, for fixed values of the 



other parameters, has the algebraically largest value of y' as defined 

R 



above. 

Incidentally, for purely sinusoidal oscillations, that is, for 



"r- s r- 0 



(4.11; 



Eq. (4.10) reduces to 



y 1 = y < 0 

1 EL r EL - 



(4.12 
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In other words the stability is the same with respect to all uniformly 

translating reference frames in this case, irrespective of the values of 

• • 

and z- Q involved. Thus the criterion of stability is simple and 
unambiguous under these circumstances. 

When we revert to the general case, however, the apparent stability 

• • 

as expressed by Eq. (4.10) depends on the velocities x^ , of the 

chosen reference axes. Thus the concept of instability is not completely 

• • 

and uniquely defined until parameters X q and z q are designated in an 
appropriate way. 

In the present context it is reasonable to require that these velo- 
cities be specified in such a way that the resulting stability criterion 
be consistent with Eqs . (4.4) and (4.5) for the two special cases to which 
these relations apply. Comparison of Eq. (4.10) with Eqs. (4.4) and (4.5) 
reveals that this will indeed be the case provided that we set 



z = 0 
o 

(4.13) 

x > 0 
o 



Consequently, Eqs. (4.9) and (4.10) now simplify to the form 



Y 



RL 



Y 



IL 



Y RL + a R X o - ° 



Y IL + a l X o 



(4.14) 



where the magnitude of X q still remains to be defined. In this connec- 
tion, a method that is followed by some investigators is to choose x 

o 

in such a way as to give 



Y 



IL 



= 0 



(4.15) 
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The corresponding value of x^ , which is now termed the phase velocity, 
is seen to be 



x 

o 




(4.16) 



and this quantity is normally positive as required. Consequently, with 
respect to this particular reference frame, the flow is stable, or at least 
neutrally stable, provided that 



a R 

v T = y — ( ) y < 0 

t RL y RL K a z J t IL - 



(4.17) 



It is also possible to define a phase velocity which satisfies Eq* 
(4.15) in a more general manner. In place of setting z q = 0 , we may re- 
quire that 

• 2 *2 

x + z = a minimum (4.18) 

o o 

However, this option contravenes Eq. (4.13) and will therefore not be con- 
sidered further. 

Despite the fact that Eq. (4.17) is the stability criterion preferred 
by some writers, this investigator does not recommend its use, for the 
following reason. Note that the restriction imposed by Eq. (4.15) is 
essentially arbitrary. It is a restriction on the angular frequency and 
as such has no necessary relation to growth or decay rate. Consequently 
the stability criterion expressed by Eq, (4. 17), while mathematically correct 
with respect to the axes to which it pertains, is as arbitrary as are 
those axes themselves. 

To eliminate the arbitrary element involved in Eq, (4.17), this investi 
gator advocates a criterion based on adopting a Lagrangian viewpoint, Con- 
sider a given streamline of the mean flow located at some specified value 
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of coordinate y . The fluid particles which lie along this streamline 
advance downstream at an average velocity U relative to fixed axes. 

We now stipulate that the flow along this streamline shall be deemed 
stable, neutrally stable or unstable, respectively, according to whether 
the particles advancing along that streamline undergo gradually in- 
creasing, unchanging or gradually decreasing amplitudes of oscillation. 
This amounts to letting the references axes move with the mean velocity 
of the fluid particles along the chosen streamline. Hence we set 



= u(y) = J (1-y 2 ) 
z = 0 



(4.19) 



Consequently, the flow along the given streamline maybe said to be 
stable, or at least neutrally stable, provided that for the least stable 
root 

Y RL = Y RL + a R - 0 (4.20) 

According to this Lagrangian criterion of stability, the least stable 

root is that root which, for given values of a , U(y) and the other 

K 

parameters, has the highest algebraic value of y' . This will also be 
the root which has the highest value of y with respect to the original 
fixed axes . 

It also follows from the foregoing concept that different streamlines 

exhibit different degrees of instability (except for the special case 

where a = 0) . Of particular importance is the least stable streamline, 

R 

that is, the streamline which, for given values of y^ and , has the 

algebraically largest value of y ' • Eq. (4,20) reveals that if a be 

RL R 
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positive, then the streamline at mid channel is the least stable, while 

if a be negative, then the streamline at the wall is the least stable. 

R 

In either case we can state that a condition of incipient instability is 
reached when 

vk m 0 (4 - 21) 



for the least stable streamline. 

From the foregoing considerations we may conclude finally that the 
flow is stable, or at least neutrally stable, 



when > 

K 



if 



Y < 

y rl - 



2 a R 



(4.22) 



when a < 0 , if y_ T < 0 

K KL 

Notice that the least stable root Yt, t in Eqs. (4.22) is now the 

KLi 

ordinary root with respect to stationary axes; the moving axes have served 
their purpose and are not needed further. Also notice that the second of 
Eqs. (4.22) is identical to the conventional stability criterion of Eq. 
(3.3). Thus only if a > 0 do we need to modify the conventional 
stability criterion in the manner shown. 

The generalized Lagrangian criterion of stability as developed in 
this section and summarized concisely in Eqs. (4.22), forms an essential 
element of the proposed research. This generalization of the stability 
criterion represents an original contribution of the present investigator 
and, so far as he is aware, does not duplicate the published work of any 
other researcher. 

Moreover, preliminary results calculated by' this method and reported 
by Harrison, Newby, and Johnston, Refs. (6), (7) and (8), appear reasonable 
and are encouraging. They show, for example, that the introduction of 
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negative values of a R is definitely destabilizing and lowers the critical 
Reynolds number below that which applies for the conventional case, 
a R = 0 . This is a significant result* Moreover, it is a result that 
coul4 not have been discovered if the investigation had been guided by 
the overly restrictive criterion of neutral stability as expressed by 
Eqs. (4.3). 

The present theory also shows that for negative values of a 

R 

instability starts at the walls at some critical Reynolds number, and with 

gradually increasing Reynolds number it spreads progressively inward. 

This is in agreement with experimental observation, Schlicting, Ref. (3). 

Application of the foregoing concepts permits us to assign the four 

parameters a , a , 3 , 3 T arbitrarily and to calculate by a process of 
R X R X 

numerical trial and correction the corresponding Reynolds number for 
incipient instability. The results disclosed by calculations of this 
kind for a suitable range of values of the above four parameters can be 
symbolized in the form 



R e ] . R e ]; *-“r’ “i’ 2 R’ 6 I-* 

Of particular interest is the question whether there exists some 
particular combination of parameters a R , a^., 3 R , 3-j. for which R^ 

iL 

attains an absolute minimum value, or whether a , a , 8,, and B T > 

R X K J- 

if otherwise unrestricted, may always be so chosen as to reduce R 



(4.23) 



below any pre-assigned limit. This amounts to asking whether, or in what 
sense, a critical Reynolds number may be said to exist below which the flow 
field is stable to all three dimensional perturbations. While this pro- 
blem has been solved for the special case where a = 3^ = 0 , it remains 
relatively unexplored for the far more complex generalized case. 
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5, Transformation of Parameters for Plane Poiseuille Flow 



The foregoing development has shown that the linearized stability 
problem of plane Poiseuille flow can be formulated in terms of the para- 
meters ct^, a^, $-£, and . In this section we introduce an 

alternative but equivalent set of parameters which are here symbolized as 
A, <$>> k and R^* . This change in parameters yields an important 
theoretical result and also permits the numerical calculations to be 
significantly reduced. 

The fundamental relations, Eqs. (2.11) and (2.12), which govern the 

solution for the eigenvalues and eigenfunctions are seen to contain two 

important complex constants as shown in the following two equations. The 

left side of each equation expresses each constant in terms of the four 

original parameters a , a , B , 3 and the right side in terms of the 

KIRI 

four new parameters A, <j> , 0 and < . Thus 



(c R + i a i ) 2 + (B r + i Bj.) 2 = A 2 e U 



(o R + 1 V 



«Ae 



( d) + Q) 



(5.1) 



Eqs. (5.1) suffice to fix the solution for the four new parameters 

A, <j>, 0 and k as functions of the four original parameters 

ot , a , g , B x or vice versa. The details of these conversions in 
R I R I 

either sense are given by Newby, Ref. (7), and are not repeated here. 
Newby also gives certain auxiliary relations associated with this trans- 
formation of parameters which need not be considered in the present 
discussion. 
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It is useful, however, to define here the important auxiliary 
parameter 

2 2 2 

\ ‘“r +6 R <5 ' 2 > 

The relationships symbolized by Eqs. (5.1) can also be summarized 

graphically in the complex plane as shown in Fig. 5.1. 

The vector OP in the diagram represents the complex quantity 
2 2 

a . The vector PQ represents 6 . The geometry of the diagram is 

completely fixed when the coordinates of the points P and Q are 

specified. Notice that the locations of these two points suffice to 

determine the four characteristics a , a , g , 8 T or A, <j>, 0 and 

KIRI 

k which fix the spatial wave form of the perturbation. 

Consider the Reynolds number which corresponds to incipient insta- 
bility for a perturbation whose spatial wave form has been specified. 
This relation is symbolized by Eq. (4.23) which we repeat here for easy 
reference, namely. 



R el R el t a R* 



a l» S R’ S I ] 



(4.23) 



A remarkable and important feature of the transformation defined by 
Eqs. (5.1) and (5.2) is that it permits the relation symbolized by Eq. 
(4.23) to be simplified significantly as explained below. 

The transformation in question makes use of a special reference 
case which, for arbitrary values of parameters A, <j> and 0 , is 
characterized by a particular value of the fourth parameter, namely, 



K K * = 1 



(5.3) 
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SPATIAL WAVE FORM PARAMETERS 



It can also be shown from Eqs. (5.1) that under these conditions the 



parameter A also takes on a corresponding reference value, namely, 
R 



A R 2 A r 2 = A 2 [|sin 0| + cos 2 (})J (5.4) 

For the above reference case, the solution of the incipient sta- 
bility problem analogous to that symbolized by Eq. (4.23) reduces to a 
form which may be symbolized as follows. 

R el * R e i* = R e i* lA > *’ 0] (5 * 5) 

The important feature of Eq. (5.5) as compared with Eq. (4.23) is 
that the number of independent parameters appearing on the right has been 
reduced from four to three, namely, to A, (j>, and 0 . Owing to the 
length and complexity of the matrix calculations involved, this reduc- 
tion from four to three independent parameters represents a very great 
reduction in the overall computational burden. 

The problem that now remains, however, is how to generalize from the 
above special reference case to the more general case in which parameters 

A, <J> and 0 remain fixed at their previous values, but in which k is 

2 

allowed to vary arbitrarily. Of course parameters A r and R also 
vary accordingly. It can be shown, however, that these variations are 
governed by the following simple analytical relations, namely. 




(5.6) 
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The detailed derivation of the relations symbolized by Eqs. (5.3) 
through (5.6) lies outside the scope of this discussion but is fully 
documented in Ref. (7). 

The general nature of the relation defined by Eq. (5.6) is displayed 

graphically in Fig. 5.2. Notice that for fixed values of A, $ and 0 , 

2 -1 
an increase in A. always implies a corresponding decrease in k 
K 

2 

and hence a decrease in R _ . Moreover, if A be allowed to increase 

el R 

without bound, then R approaches zero in the limit. On the other 
2 

hand if A is limited to some specified upper bound, then R cannot 
K el 

drop below some corresponding lower bound. 

The foregoing discussion is based on holding parameters A, cf) and 0 

2 

fixed and considering the effects of possible changes in A . Let us 

R 

2 

next consider the case where A is held fixed and parameters A, <f> 

R 

and 0 are free to vary. It follows from Eqs. (5.5) and (5.6) that 

there now exists some definite combination of values for A, cf> and 0 , 

respectively, for which R^ attains an absolute minimum value. We 

term this the critical Reynolds number R ) . It is evident that 

e cr 

R^) is some definite function of the single remaining independent 
2 

parameter, A^ . One important and specific objective of the proposed 
research is to compute this function . 

2 

Eqs. (5.5) and (5.6) show that for an assigned value of A , the 

R 

critical combination of values of A, <j> and 0 which minimize the value 
of cannot be found analytically, but must be determined by a syste- 

matic process of numerical exploration over the four dimensional domain 

* 

whose coordinates are A, <£, 0 and R^ « Some preliminary numerical 
calculation of this kind has already been accomplished on a limited scale, 
and research support is needed for the purpose of extending this exploration 
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** soo -i^r) 
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I 



in a systematic way. This task forms an essential part of determining 

o 

the critical Reynolds number as a definite numerical function of A 

R 

The theory developed in this section shows that, even before any 

specific numerical calculations are undertaken, the so-called critical 

Reynolds number is not some absolute constant as it is commonly conceived 

2 

to be, but rather some definite function of A . This is a significant 

R 

new concept and represents an original contribution by the present 
investigator. 

For the special case where 

* R 2 = 0 (5.7) 

the present theory and transformations reduce to the conventional and well 
known relations associated with Squire f s theorem. This investigator knows 
of no published result, however, which is fully equivalent to the generalized 
form of these transformation as outlined in this section. These results 
are believed to represent an original contribution by the present investigator. 

The analysis of the generalized criterion of stability as presented in 
the previous section has shown that the conventional constraint expressed by 
Eq. (5.7), which is widely regarded as a necessary condition for neutral sta- 
bility, is not in fact necessary at all. The analysis of this section now 

2 

shows that positive values of A have the further effect of decreasing 

R 

the critical Reynolds number. This result helps explain why the conventional 
analysis tends to predict values of critical Reynolds number which are some- 
what too high and how slight differences in experimental set-up or test 

2 

conditions, which might be associated with different values of A , could 

h 

produce apparent inconsistencies in the observed critical Reynolds number* 

A significant point is that this explanation of discrepancies is still 



5-7 



based on a relatively straightforward linearized analysis whereas some 
investigators, failing to find any simple explanation of this kind, have 
felt compelled to resort to non-linear analyses, with all of the enormous 
complexities they entail, in an effort to make sense out of the experi- 
mental observations. 
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6. Stability of Pipe Flow 

While the stability of plane Poiseuille flow as discussed in the 
previous sections of this report is of fundamental interest because of the 
basic simplicity of the three dimensional rectilinear geometry involved, 
this configuration is of somewhat limited practical importance, and it 
does not lend itself readily to experimental study. In fact the amount of 
experimental data available for this case is extremely limited. 

Pipe flow on the other hand, while involving a geometry which is only 
slightly more complex, is of the greatest practical importance. Also, 
there exists an extensive body of experimental data on the stability of 
pipe flow. Moreover, the discrepancy between the conventional theory and 
experimental observation is extreme in this case in that the theory fails 
completely to predict the experimentally observed instability. These are 
strong reasons for generalizing the theoretical analysis for pipe flow to 
include three dimensional effects in a manner analogous to that previously 
discussed in connection with plane Poiseuille flow. Preliminary calcula- 
tions based on such a generalized analysis, while of very limited scope , 
have succeeded for the first time in disclosing theoretical instabilities 
in pipe flow . This is a noteworthy accomplishment and warrants a sustained 
follow up. 

Broadly speaking, the pipe flow problem can be analyzed by much the 
same method previously explained in connection with plane Poiseuille flow. 
Only certain details are different in the two cases. An extensive analysis 
for pipe flow is given by Harrison, Ref. (6) and Johnston, Ref. (8), and 
is not repeated here. However, we do point out in this section the main 
respects in which the analysis of the stability of pipe flow differs from 
that of plane Poiseuille flow. 
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All equations are nondimensionalized by utilizing as dimensional 

reference parameters the density p of the fluid, the radius a of the 

pipe and the mean volumetric velocity U of the fluid. 

Cylindrical coordinates x, r, 9 are used with x being in the 

direction of the flow, r being in the radial direction and 0 being the 

angular coordinate. Dimensionless coordinate r varies from zero at the 

axis to +1 at the pipe wall. Symbols e , e , e denote unit vectors 

x r 0 

in the three coordinate directions. 

It is well known that the dimensionless velocity U of the mean flow 
along the x axis is given by the simple expression 



U = 2(1 - r ) 



( 6 . 1 ) 



It is of interest to compare this relation for pipe flow with the corres- 
ponding relation for plane flow as given by Eq. (2.1). Notice that the 
maximum dimensionless velocity now equals 2 as compared with 3/2 for plane 
flow. 

The vector potential of the perturbation is chosen in the form 



W = [e (o) + e G (r) 
n x r n 



+ e 0 H n (r)] 



( 6 . 2 ) 



where the complex constants 



a and y may be expanded further as follows, 



a = a + i a 
R I 



Y = Y r + i Yj 



(6.3) 



Notice that the coefficient of 0 in the exponent in now taken not 
as a general complex quantity g , but rather as the purely imaginary 
quantities i n where n is an integer. Thus 
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(6.4) 



6 - 3^ + i Bj - 0 + i n 

6 r = 0 

= n n = 0, 1, 2, 3, ... 

\ 

The reasons for restricting the complex coefficient g in this way is 
that the velocity components, and hence the vector potential itself, must 
be strictly periodic with respect to angle 0 . This is a fortunate circum- 
stance since it eliminates parameter g^ and reduces g^ to an integer, 
thus simplifying the problem. The spatial form of the perturbation is 

therefore completely defined by the three parameters a , a and n . 

R 1 

As Refs, (6) and (8) show in full detail, the vorticity transport 
equations that result from the use of the above vector potential are found 
to reduce to the following form 



'D 4 G 



n 



, 3 
D G 



n\ 



'D 2 G 



n 



DG 



n 



)d 4 h ) + f M 3'] | D 3 h f + f M 2 '] |d 2 H ) + (DH 

n 7 vn y ^n 



+ [M o'l IB 



n 



= T 



n , 



/ (d 2 g ) (dg ) 



'] 



n 



IH 



(6.5) 



n 



where symbol D now denotes the differential operator d/dr . 

The eight matrices [M’] and [N'] which appear in Eq. (6.5) each have 

two rows and two columns, or 32 elements in all, of which 10 happen to be 

zeros. All elements are defined in detail in Refs. (6) and (8). These 

definitions are too lengthy to be repeated here<, It suffices to remark 

that all of these matrix elements are known complex functions of the three 

spatial wave form parameters a , a and n , of Reynolds number R , 

R 1 e 
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and of the radial coordinate r . Moreover, the coefficients of the 

3 4 

quantities D and D in Eqs. (6.5) are zeros, so that these equation 

are of only second order in G whereas they are of fourth order in H 

n n 

Thus only two boundary conditions are needed which involve G^ or its 

derivatives, while four are needed which involve H or its derivatives. 

n 

An interesting feature of Eqs. (6.5) is that, in general, these two 

equations are coupled, that is, both functions G and H and their 

n n 

various derivatives appear in both equations. Hence the equations must 
be solved simultaneously. For the special case n = 0 , however, these 
equations uncouple. One of them becomes an eigenvalue equation of fourth 
order in the single unknown function H q . The other becomes an independent 
eigenvalue equation of second order in the single unknown function G^ . T1 
the solution for n = 0 breaks down into two independent families of eigen- 
values and eigenfunctions. It can be shown that one of these families 
represents perturbations in the meridional plane 0 = constant , the other 
represents perturbations in the cross-sectional plane x = constant . On t 1 

other hand for n = 1, 2, 3 ... the eigenfunctions G and H which 

n n 

correspond to a particular eigenvalue remain coupled and therefore represent 
perturbations of a truly three dimensional character. 

Solution of the basic vorticity transport relations, Eqs. (6.5) must 
be found subject to the appropriate boundary conditions at the wall., r = 1 

I 

and at the axis, r = 0 . The boundary conditions are discussed in some 

detail in the next section of this report. 

If numerical values of ol., a T , n and R be specified, the real 

K i e 

part of the least stable root can be found by solution of the vorticit 

transport equations. The stability can then be determined by a Lagrangian 
criterion analogous to Eq. (4.22). In other words the flow is stable, or 
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at least neutrally stable. 



when 


“r* 0 • 


..if 


% - ' 2 “r 


when 


“r - 0 ’ 


if 


YrIO 



( 6 . 6 ) 



Applications of the foregoing concepts permits us to assign the three 
parameters a R , and n arbitrarily and to calculate by a process of 

numerical trial and correction the corresponding Reynolds number for 
incipient instability. The results disclosed by calculation of this kind 
for a suitable range of values of the above three parameters can be 
symbolized in the form 



R el " R el >V V nl 



(6.7) 



Of particular interest is the question whether, if we hold a con- 

R 

stant at some arbitrary value, there exists some particular combination of 

parameters a^. and n for which R^ attains an absolute minimum value* 

We term this the critical Reynolds number R ) . It is evident that 

e cr 

R^) some definite function of the single remaining independent para- 

meter, cip . One important and specific objective of the proposed research 
is to compute this function . 

While the stability of pipe flow has been extensively investigated for 
the special case a = 0 , it nevertheless remains relatively unexplored 
for the more general case in which a ^ 0 . Moreover, in pipe flow certain 
boundary conditions at the axis become so problematical as to raise serious 
questions about the validity of any purported solution that does not properly 
take these subtleties into account* These aspects are considered in the 
next two sections* 
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7. Boundary Conditions for Pipe Flow- 

In pipe flow, the final vorticity transport equations are of second 

order in function G and of fourth order in function H . Hence two 

n n 

boundary conditions are needed with respect to G and four with respect 

n 

to H It turns out that three of these conditions are defined at the 



n 



wall, r = 1 , and the remaining three at the axis, 
The boundary conditions at the wall, r = 1 , 
no slip restriction. Thus, for all values of index 

u n (l) -O' 

v (1) = 0 ; 
n 



r = 0 . 

are based on the usual 
n 



(7,1) 



w (1) = 0 
n 



Inasmuch as the governing differential equations are finally expressed 

in terms of the two basic functions G (r) and H (r) , it is useful to 

' n n 

rewrite Eqs. (7,1) in like terms. In this connection recall that the 
perturbation velocity v^ is related to the vector potential as 

follows. 



v = VxW = 
n n 



-+ 

e u (r) 
x n 



+ e v (r) 
r n 



e e V r) 



(ax + in© + yt) 



(7.2) 



From Eqs. (6.2) and (7.2) we readily find that the three velocity com- 
ponents are given by the following expressions. 



u 

n 



= DH 

n 





in G ) 
n 



v = ~ aH 
n n 



w = 4* aG 
n n 



(7.3) 



7-1 



It is also useful in the present context to expand functions 

and H in Maclaurin series. Thus 
n 



G 

n 



- V o) 



DG (o) 
n 



r + 



d 2 g 

n 



(•>£ 



t f I 



H 



n 



■ V o) 



DH (o) 
n 



r + 



D 2 H (o) 
n 



2 

r 

21 



G 

n 



(7.4) 



With the aid of Eqs . (7.3) and (7.4) it becomes a straightforward 
task to rewrite Eqs. (7.1) in the required form. The following results 
are obtained which apply for all values of index n , namely, 

G n (l) = 0 

H n (l) - 0 (7.5) 

DH (1) = 0 
n 

Derivation of the remaining three boundary conditions which apply at 

the axis r = 0 is considerably more complicated. We start by observing 

that Eq. (7.2) correctly gives the perturbation velocity v^ at any point 

which is not on the axis but does not necessarily define a meaningful 

velocity at a point on the axis. To express the velocity on the axis it is 

not sufficient merely to set r = 0 in Eq. (7.3). The reason is that 

-*• in© 

Eq. (7.2) contains the quantities e^, eg and e which are functions 
of © , whereas coordinate 0 , while determinate for any point off the 
axis, is indeterminate at the axis itself. 

We may circumvent this difficulty by starting at some arbitrary point 
near but not on the axis, and proceeding toward the axis along a radial line 
for which 0 remains constant. In the limit as r 0 , the final location 
of the point on the axis is independent of the coordinate 0 of the line 
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along which that point is approached. Likewise, the velocity v n at that 
point is independent of the coordinate 0 which characterizes the path of 
approach. Hence we may write 



In applying 
is a constant but 



9 / lim 

30 (r+0 



Eq. (7.6) to 

e and e n 
r 0 



v n (x, r, 0, t)| = 0 (7.6) 

Eq. (7.2) we make use of the fact that e 

x 

are functions of 0 , so that 




= 0 





(7.7) 




e 

r 



It now 
arbitrary 0 



follows easily from Eqs . (7.2), (7.6) and (7.7) that for 
and arbitrary n 



n u (o) = 0 



i n v (o) - w (o) =0 (7.8) 

n n 



v (o) + in w (o) = 0 
n n 



Eqs. (7.8) stipulate the conditions that must be satisfied in order 
that Eq. (7.2) shall define a unique and determinate value of in the 

limit as r 0 « The value of 0 at this limit may be chosen arbitrarily 
without affecting the resulting value of v q . 

Eqs. (7,8) break down into three distinct cases, namely, n = 0 , 
n = 1 and n > 2 . Thus 
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For n = 0 



v (o) = 0 
o 

w (o) = 0 
o 

Notice that u (o) does not vanish in this case, 
o 

For n = 1 

u 1 (o) = 0 

v 1 (o) + i w^o) = 0 



(7.9) 



(7.10) 



Notice that v^(o) and w^(o) 
not vanish in this case. 

For n=2, 3, 4, 5, ... 

u (o) = 0 
n 

v (o) = 0 
n 

w (o) = 0 
n 



although mutually constrained, do 



(7.11) 



It is significant that u 0 (°) > v (o) and w^(o) are the only non- 
zero velocity components at the axis. Thus the cases n = 0 and n = 1 
suffice to account fully for the three non zero perturbation velocity com- 
ponents that can occur at this location. 

It is now useful to express Eqs . (7.9), (7,10) and (7,11) in terms of 

functions G and H . The following results are obtained 
n n 



For n = 0 



G (o) = 0 
o 



H (o) = 0 
o 



(7,12) 



Note that one additional boundary condition is still needed for this 



case , 
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For n = 1 



H^o) - i G 1 (o) = 0 
2DH (o) - i DG (o) = 0 



(7.13) 



Note that one additional boundary condition is still needed for this 
case also. 

Notice in addition that this case is exceptional in that it represents 

the only value of n for which H (o) and G (o) do not vanish, 

n n 

For n = 2 , 3, 4, 5, 



G (o) = o 
n 

H n (o) = 0 (7,14) 

2DH (o) - i n DG (o) = 0 
n n 

It is now clear that for n ^ 2 , Eqs. (7.5) and (7.14) fix all six 
of the boundary conditions needed to define a determinate solution. On the 
other hand the cases n = o and n = 1 each still lack one final necessary 
boundary condition. 

The one remaining boundary condition needed for the case n = 0 can 
readily be deduced from well known physical laws, Consider a small fluid 
element of cylindrical shape, of length 6 x and radius Sr , whose axis 
coincides with the x coordinate axis. Let denote the longitudinal 

shear stress acting on the curved surface of the element, let p be the 
pressure, p the density and (— ) the axial acceleration. We may write 
the equation of motion in the x direction for this element as follows. 



T 27 r 6 r 5x - (|^) Sxtt (Sr ) 2 = p 7 r(Sr ) 2 5 x (-— 7 ) 
xr ox at 



(7.15) 
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Proceeding to the limit 6r ** 0 , we observe that the last two terms 



of Eq. (7.15) are of higher order than the first. Hence 



lim x = lim t 



= 0 



(7.16) 



xr 



xr 



6r ^ 0 r 0 



It then follows from the known relation between stress and strain 
rate for a viscous fluid that the corresponding shear strain rate also 
vanishes along the axis. It can readily be shown that this requires in 
turn that 



making use of Eqs. (7.9), we find that the vanishing of the shear strain 
rate requires finally that 



This fixes the last boundary condition needed for the case n = 0 . 
All of the boundary conditions required for the pipe flow problem 
are therefore determined above except one. This is the last boundary con- 
dition needed to define a determinate solution for the important case 
n = 1 . 

Unfortunately, this last required boundary condition is also the most 

difficult to formulate in a clear, unequivocal and rigorous manner. The 

axis represents a singular line along which many of the usual analytical 

expressions for quantities like v, ca, Vxw and so on break down. The 

functions G and H and their various derivatives, especially the 
n n 

higher derivatives, are not necessarily continuous at r = 0 , Conse- 
quently, intuitive notions based on various symmetry relations and 




(7.17) 



Upon expressing this in terms of functions G q and and upon 



D 2 H (o) = 0 

O 



(7.18) 
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continuity assumptions become unreliable guides. 

We call attention at this point to the analogy that exists between 
perturbation vorticity and perturbation velocity v in this flow field. 

Thus 

. (ax + in© + yt) 

“n = V X ^n = l®x + 8 r n ^ r) + 8 0 8 (7.19) 



The analogy between this relation and Eq. (7.2) implies that if 
Eq. (7.19) is to remain valid at r = 0 , then we must have, by analogy 
with Eqs . (7.8), 



n £ n (o) = 0 

i n n(o) - c (o) = 0 
n n 

n (o) + in C (o) = 0 
n n 

In particular for n = 1 , these reduce to 

? L (o) = 0 
n 1 (o) + i ^(o) = 0 



(7.20) 



( 7 . 21 ) 



Of course Eqs. (7.13) already supply two of the three needed boundary 
conditions, and Eqs. (7.21) are under consideration here only in connection 
with the one additional constraint that is still lacking. If we invoke the 
first of Eqs. (7.21) for this purpose, the following constraint is obtained, 
namely. 



DH 1 (o) - i 2DG 1 (o) = 0 



(7.22) 
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On the other hand, if we invoke the second of Eqs. (7.21) for the same 
purpose, a quite different result is obtained, namely 

D 2 H 1 (o) - i D 2 G 1 (o) = 0 (7.23) 

Unfortunately, with only one degree of freedom still at our disposal, 
there is no way to satisfy both Eq. (7.22) and Eq. (7.23) simultaneously. 

Hence both of Eqs. (7.21) cannot be satisfied simultaneously either. This 
simply means that for the case in question, n = 1 , Eq. (7.19) inevitably 

| as 

breaks down at r = 0 . Moreover, the attempt to utilize the analogy between 

I 

Eq. (7.19) and (7.2) leads to an overspecification of the needed final 
boundary condition . 

Further study is needed in order to reformulate the required final 
boundary condition at the axis in a way that is unambiguous and free of 
arbitrary assumptions. This is an important question, and the difficulties 
involved in correctly formulating this particular boundary condition might 
perhaps account for the failure of the conventional theory to disclose the 
known instability of pipe flow. A cricial review of this question will form 
an important part of the proposed research . 

For the present, we shall tentatively assume that the first of Eqs. (7.21 
provides the appropriate final boundary condition, and that the second of 
these equations may be discarded. This amounts to retaining Eq. (7.22) and 
discarding Eq. (7.23). 

Upon combining Eqs. (7.13) and (7.22), we find that the three final 
boundary conditions at the axis simplify for this case to the following form. 
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For n = 1 



DG (o) = 0 

H 1 (o) - i G 1 (o) = 0 (7.24) 

DH^o) = 0 

The final overall boundary conditions for all values of n are now 
summarized in Eqs. (7.25) through (7.28) of Table 7.1. 
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Table 7,1 Boundary Conditions fox Pipe Flow 
For n = 1 



G o (l) = 0 

(7.25) 


H o (l) = 0 


G (o) = 0 
o 


DH (1) = 0 

(7.26) 

H (o) = 0 
0 




D 2 H (o) = 0 
o 



This case is exceptional in that the solutions for functions G^ 



and H become uncoupled for n = 

n 

For n = 1 


= 0 . 


0 
II 

N 

1 1 

r— 1 

o 


o 

o 

II 

o 


H^l) = o 


H-^o) - i G 1 (o) = 0 (7.27) 


DH (1) = 0 


DH 1 (o ) = 0 



This case is exceptional in that the quantities G (o) and H (o) 

n n 



become non-zero for n = 1 . 

For n=2, 3, 4, 5, .... 




G n (l) - 0 


G (o) = 0 
n 


H n (l) - 0 


H q (o) = 0 (7.28) 


DH (1) = 0 
n 


2DH (o) - in DG (o) = 0 
n n 
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3. Radial Symmetry Relations and Continuity Conditions at Axis 

It is quite common in dealing with a variety of physical problems 
vhich involve a cylindrical geometry to derive certain boundary conditions 
at the axis from considerations of symmetry. Hence it is natural to in- 
quire whether any such constraints apply also to the pipe flow problem. 

The main purpose of the discussion in this section is to show that although 
the solution of the present problem does indeed satisfy certain relations 
of radial symmetry, the attempt to derive the boundary conditions on this 
basis involves certain hidden pitfalls. We may even conjecture that such 
pitfalls might be a factor in the failure of conventional efforts to find 
any theoretical instability in pipe flow, although we know from experience 
that such instability does in fact occur. 

With reference to symmetry, consider an arbitrary point P with 
coordinates x, r, 0 . A diametrically opposite point P f may now be 
located with coordinates x, r, (0 + it) . While negative values of coordi- 
nate r are not ordinarily used, we can elect in this instance to employ 
them, and to designate the coordinates of point P* in the alternative 
but equivalent form x, -r, 0 . Of course the perturbation velocity at 
point P* , call it v , must be the same regardless of which of the above 
two alternative coordinate descriptions are employed. Thus 

v [x, -r, 0, t] = v [x, +r, (0 + tt) , t ] (8.1) 

Harrison, Ref. (6), has shown that as a consequence of the relation 
symbolized by Eq. (8.1), the functions and H^ must satisfy symmetry 

constraints of the form 
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( 8 . 2 ) 



G (™r) = - (-l) n G (+r) 
n n 

H (-r) = - (-l) n H (+r) 
ti n 

This tells us that if n be even, functions G and H have odd 

n n 

symmetry with respect to r , while if n be odd, G and H have even 

n n 

symmetry. 

Consequently, for n even, 

G (-r) = - G (+r) 
n n 

D 2 G (-r) = - D 2 G (+r) (8.3) 

n n 

D 4 G (-r) = - D 4 G (+r) etc. 
n n 



and similarly for function H and its even derivatives. 

n 

Likewise, for n odd, 

DG (-r) = - DG (+r ) 
n n 

D 3 G (-r) = - D 3 G (+r) etc. 
n n 



(8.4) 



and similarly for the odd derivatives of function H 

n 

If we now assume that the functions and derivatives in question are 
also continuous at r = 0 , then the above relations have the following 
further consequences.* 



For n even 

G (o) = 0 
n 

D 2 G (o) = 0 
n 

D 4 G (o) = 0 
n 

etc. 



H (o) = 0 
n 

D 2 H (o) = 0 (8.5) 

n 

D 4 H (o) = 0 
n 

etc. 
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For n odd 



DG (o) = 0 



n 



DH (o) = 0 



n 



( 8 . 6 ) 




n 




etc. 



etc. 



An error that is sometimes made is to assume that the known symmetry 
relations expressed by Eqs. (8.2) necessarily imply the consequences 
summarized in Eqs. (8.5) and (8.6). In fact, however, these consequences 
apply if and only if the quantities in question are known to be continuous 
at r = 0 . Inasmuch as the mathematical functions ordinarily encountered 
in connection with incompressible flows are usually highly continuous, it 
is easy to overlook the fact that Eqs. (8.5) and (8.6) actually involve a 
hidden assumption of continuity that may or may not be justified in the 
present instance. 

The situation is further clouded by the fact that a number of the 
boundary conditions derived in the previous section from entirely different 
considerations, and summarized in Table 7.1, happen to coincide with the 
symmetry/continuity relations of Eqs. (8.5) and (8.6). These are the 
following 

For n = 0 



G (o) = 0 
o 



H (o) = 0 
o 

D 2 H (o) = 0 
o 



(8.7) 



For n = 1 



DG (o) = 0 
DH (o) = 0 



( 8 . 8 ) 



For n = 2, 4, 6, , . , 

G (o) = 0 
n 

H (o) = 0 
n 



(8.9) 



Moreover, for the case n = 3, 5, 7, . . . , if Eqs. (8,6) could be 
assumed to apply, they would also ensure that the condition 

2DH (o) - in DG (o) = 0 (8,10) 

n n 

is satisfied identically. 

The correct interpretation of Eqs. (8.7), (8.8) and (8.9) would seem 
to be that these relations, derived earlier on quite independent grounds, 
ensure that the corresponding particular quantities must also be continuous 
at r = 0 . But this restricted conclusion does not warrant the indiscri- 
minate extension of Eqs. (8.5) and (8.6) to other cases without further 
justification. 

Furthermore, there are a number of boundary conditions listed in 
Table 7.1 which, while not necessarily incompatible with Eqs, (8.5) and 
(8.6), are entirely unrelated to the latter. These are the following 

For n = 1 

H 1 (o) - i G 1 (o) - 0 (8,10) 

For n = 2, 4, 6, ... 

2DH (o) - i n DG (o) = 0 (8.11) 

n n 
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For n = 1, 3, 5, , , . 

G (o) = 0 

( 8 . 12 ) 

H n (o) = 0 

If we proposed for some reason to derive the boundary conditions 
exclusively from symmetry relations and continuity assumptions, it would 
be necessary to abandon the conditions listed in Eqs. (8.10), (8,11) and 
(8.12), and to substitute for them other relations drawn from Eqs. (8.5) 
(8.6). At the present time no adequate physical or mathematical justifi- 
cation can be offered for such a change. 

In an effort to clarify some of the conceptual complications 
associated with Eqs, (8,5) and (8,6), Johnston in Ref. (8) investigated 
the constraints that would theoretically have to be satisfied in order to 
permit the direct substitution r = 0 into the vorticity transport equa- 
tions, His principal result pertains to the case n >_ 2 . He found for 

this case that functions G and H and their first four derivatives 

n n 

would all have to vanish at r = 0 to allow this. His result involves 
ten constraints in all. Unfortunately, with only three degrees of freedom 
available at the axis, only three of these ten constraints can actually be 
enforced. While the details are different for n = o and n = 1 , the 
number of constraints that would theoretically be required for this pur- 
pose always exceeds the three available degrees of freedom. Hence we 
conclude that the simple substitution r = 0 into the vorticity transport 
equations, in their original analytical form, is never permissible. In 
other words, the axis remains irreducibly a singular point, and Johnston’s 
analysis, while correct in itself, does not provide the three and only 
three boundary conditions needed at that point. 
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In view of the various considerations discussed above, it is 
believed that Table 7.1 provides the best overall formulation of the 
boundary conditions that is available at this time. 
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9. Status of Pipe Flow Stability Calculations 



Harrison, Ref. (6), correctly formulated the basic equations for 
pipe flow and the boundary conditions for n = 0 and n = 1 . Unfortu- 
nately, he was unable to compute results for these cases owing to a minor 
error in his computer program. Johnston, Ref. (8), corrected this error 
and computed valid results for n = 0 , but the scope of his calculations 
was very limited. Johnston, in turn, formulated the boundary conditions 
incorrectly for all cases other than n = 0 . Hence a specific aim of 
the proposed research is to amend and extend the above calculations, 
using the corrected boundary conditions of Table 7.1. It is felt that 
the revised boundary conditions for n = 1 are particularly likely to 
disclose instabilities that have not appeared before. 

Despite the error in boundary conditions involved for some of the 

cases, Johnston’s results in Ref. (8) suffice to indicate the strong 

probability that small negative values of a are definitely destabilizing. 

An important aim of the proposed research is to explore this question more 

fully and adequately. An even more exciting possibility is that with the 

revised boundary conditions discussed above for n = 1 , instability might 

possibly be found even for a = 0 . If so, this result would represent 
R 

a very significant achievement . 
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10. Summary of Research Objectives 



The first general objective of the research outlined in this report 
is to carry out the calculations described at least to the extent necessary 
to provide significant support for the overall validity of the present 
theory and for the various innovations that it embodies. For the case of 
plane flow this entails exploring some typical portion of the boundary of 
incipient instability as expressed symbolically in the form 

R ei* = R ei* [A ’ 01 (10.1) 

For the case of pipe flow, it entails a corresponding exploration of the 
form 

R el = R el ^ a R’ a I» (10.2) 

A particular objective of the work on pipe flow is to confirm the 
adequacy of the revised boundary conditions at the axis as summarized in 
Table 7.1. The case n = 1 warrants special attention as having possibly 
the greatest potential for disclosing instabilities not detected previously. 

The ultimate objective is to carry the above calculations forward to 
the extent that it becomes possible to plot the two important functions 
symbolized below for plane flow and pipe flow, respectively, namely. 



R = R (X„) plane flow 

e e R 

cr cr 

R = R (a R ) pipe flow 

cr cr 



(10.3) 

(10.4) 



The demonstration of the very existance of these two functions would 
represent a noteworthy and fundamental new contribution to the theory of 
hydrodynamic stability. 
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An essential collateral objective at all stages of the research is 
to report progress periodically in the technical literature. Enough pro- 1. 

gress has already been achieved to warrant at least one initial technical 
paper of this kind. Support is needed as much for preparing such a paper 
on progress already made as for continuing further analysis and calculations. 1 2 
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